Draft version February 5, 2008 

Preprint typeset using UNgK. style emulateapj v. 1 1/12/01 



A BLIND ESTIMATION OF THE ANGULAR POWER SPECTRUM OF CMB ANISOTROPY FROM WMAP 

Rajib Saha 12 , Pankaj Jain 1 , Tarun Souradeep 2 

Draft version February 5, 2008 

ABSTRACT 

Accurate measurements of angular power spectrum of Cosmic Microwave Background (CMB) radiation has 
lead to marked improvement in the estimates of different cosmological parameters. This has required removal of 
foreground contamination as well as detector noise bias with reliability and precision. We present the estimation of 
CMB angular power spectrum from the multi-frequency observations of WMAP using a novel model-independent 
method. The primary product of WMAP are the observations of CMB in 10 independent difference assemblies 
(DA) that have uncorrelated noise. Our method utilizes maximum information available within WMAP data by 
linearly combining all the DA maps to remove foregrounds and estimating the power spectrum from cross power 
spectra of clean maps with independent noise. We compute 24 cross power spectra which are the basis of the 
final power spectrum. The binned average power matches with WMAP team's published power spectrum closely. 
A small systematic difference at large multipoles is accounted for by the correction for the expected residual 
power from unresolved point sources. The correction is small and significantly tempered. Previous estimates have 
depended on foreground templates built using extraneous observational input. This is theirs? demonstration that 
the CMB angular spectrum can be reliably estimated with precision from a self contained analysis of the WMAP 
data. 

Subject headings: cosmic microwave background - cosmology: observations 



1. INTRODUCTION 

Remarkable progress in cosmology has been made due 
to the measurements of the anisotropy in the cosmic mi- 
crowave background (CMB) over the past decade. The ex- 
traction of the angular power spectrum of the CMB anisotropy 
is complicated by foreground emission within our galaxy 
and extragalactic radio sources, as well, as the detector 

! noise (Bouchet & Gispert 1999; Tegmark & Efstathiou 1996). 

■ It is established that the CMB follows a blackbody distri- 
bution to high accuracy, (Mather et al. 1994 & 1999). Hence, 
foreground emissions may be removed by exploiting the 
fact that their contributions in different spectral bands are 
considerably different while the CMB power spectrum is 
same in all the bands (Dodelson 1997; Tegmarkef a/. 2000a; 
Bennett et al. 2003b; Tegmark 1998). Different approaches 
to foreground removal have been proposed in the liter- 

' ature (Bouchet & Gispert 1999; Tegmark & Efstathiou 1996; 
Hobsonefa/. 1998) 

(Maino et al. 2002,2003; Eriksen et al. 2005) 

The Wilkinson Microwave Anisotropy Probe (WMAP) ob- 
serves in 5 frequency bands at 23 GHz (K), 33 GHz (Ka), 
41 GHz (Q), 61 GHz (V) and 94 GHz (W). In the first 
data release, the WMAP team removed the galactic fore- 
ground signal using a template fitting method based on a 
model of synchrotron, free free and dust emission in our 
galaxy (Bennett et al. 2003b). The sky map around the galac- 
tic plane and around known extragalactic point sources were 
masked out and the CMB power spectrum was then obtained 
from cross power spectra of independent difference assem- 
blies in the 41 GHz, 61 GHz and 94 GHz foreground cleaned 
maps (Hinshaw et al. 2003a). 

A model independent removal of foregrounds has been pro- 
posed in the literature 

(Tegmark & Efstathiou 1996). The method has also been im- 
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plemented on the WMAP data in order to create a foreground 
cleaned map (Tegmark et al. 2003). The main advantage of 
this method is that it does not make any additional assump- 
tions regarding the nature of the foregrounds. Furthermore, the 
procedure is computationally fast. The foreground emissions 
are removed by combining the five different WMAP bands by 
weights which depend both on the angular scale and on the lo- 
cation in the sky (divided into regions based on 'cleanliness'). 
However, this analysis did not attempt to remove the detec- 
tor noise bias (Tegmark et al. 2003). Consequently, the power 
spectrum recovered from the foreground cleaned map has a lot 
of excess power at large multipole moments due to amplifica- 
tion of detector noise bias beyond the beam resolution. 

The prime objective of our paper is to remove de- 
tector noise bias exploiting the fact that it is uncor- 
related among the different Difference Assemblies (DA) 
(Hinshaw et al. 2003a; Jarosik et al. 2003). The WMAP 
data uses 10 DA's (Bennett et al. 2003a; Bennett et al. 2003c; 
LimonefaZ.; Hinshaw et a/.2003b ), one each for K and Ka 
bands, two for Q band, two for V band and four for W band. 
We label these as K, Ka, Ql, Q2,V1, V2, Wl, W2, W3, and 
W4 respectively. We eliminate the detector noise bias using 
cross power spectra and provide a model independent extrac- 
tion of CMB power spectrum from WMAP first year data. So 
far, only the 3 highest frequency channels observed by WMAP 
have been used to extract CMB power spectrum and the fore- 
ground removal has used foreground templates based on extrap- 
olated flux from measurements at frequencies far removed from 
observational frequencies of WMAP (Hinshaw et al. 2003a; 
Fosalba & Szapudi 2004; Patanchon et al. 2004). We present a 
more general procedure where we use observations from all the 
5 frequency channels of WMAP and do not use any extraneous 
observational input. 

2. METHODOLOGY 
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2. 1 . Foreground Cleaning 

Up to the foreground cleaning stage, our method is similar 
to Tegmark & Efstathiou (1996) and Tegmark et al. (2003). In 
Tegmark et al. (2003), a foreground cleaned map is obtained 
by linearly combining 5 maps corresponding to one each for 
the different WMAP frequency channels. For the Q, V and 
W frequency channels, where more than one maps were avail- 
able, an averaged map was used. However, averaging over the 
DA maps in a given frequency channel precludes any possi- 
bility of removing detector noise bias using cross correlation. 
In our method we linearly combine maps corresponding to a 
set of 4 DA maps at different frequencies. We treat K and Ka 
maps effectively as the observation of CMB in two different 
DA. Therefore we use K and Ka maps in separate combina- 
tions. In case of W band 4 DA maps are available. We simply 
form an averaged map taking two of them at a time and form 
effectively 6 DA maps. Wij represents simply an averaged map 
obtained from the z* and f h DA of W band. (Other variations 
are possible. We defer a discussion to a more detailed publi- 
cation (Saha et al. 2006)). In table 1 we list all the 48 possible 
linear combinations of the DA maps that lead to 'cleaned' maps, 
Ci and CAi's, where i = 1, 2, . . . , 24. 



(K,KA)+Q1+V1+W12= 
(K,KA)+Q1+V1+W13= 
(K,KA)+Q1+V1+W14= 
(K,KA)+Q1+V1+W23= 
(K,KA)+Q1+V1+W24= 
(K,KA)+Q1+V1+W34= 
(K,KA)+Q2+V2+W12= 
(K,KA)+Q2+V2+W13= 
(K,KA)+Q2+V2+W14= 
(K,KA)+Q2+V2+W23= 
(K,KA)+Q2+V2+W24= 
(K,KA)+Q2+V2+W34= 



<C1,CA1) 

<C2,CA2) 

<C3,CA3) 

<C4,CA4) 

<C5,CA5) 

<C6,CA6) 

<C7,CA7) 

<C8,CA8) 

<C9,CA9) 

<C10,CA10) 

<C11,CA11) 

(C12.CA12) 



(K,KA)+Q1+V2+W12= 
(K,KA)+Q1+V2+W13= 
(K,KA)+Q1+V2+W14= 
(K.KA)+Q1+V2+W23= 
(K.KA)+Q1+V2+W24= 
(K,KA)+Q1+V2+W34= 
(K,KA)+Q2+V1+W12= 
(K.KA)+Q2+V1+W13= 
(K,KA)+Q2+V1+W14= 
(K,KA)+Q2+V1+W23= 
(K.KA)+Q2+V1+W24= 
(K.KA)+Q2+V1+W34= 



(C13.CA13) 
(C14.CA14) 
(C15.CA15) 
(C16.CA16) 
(C17.CA17) 
(C18.CA18) 
(C19.CA19) 
(C20.CA20) 
(C21.CA21) 
(C22.CA22) 
(C23.CA23) 
(C24.CA24) 



Table 1 

Different combinations of the DA maps used to obtain 

THE FINAL 48 CLEANED MAPS. 

Following the approach of Tegmark & Efstathiou (1996), we 
introduce a set of weights, [Wi] = (w] ,w 2 ,w ; 3 ,wf) for each of 
the 4 DA in the combination, which defines our cleaned map as 
the linear combination 



Clean 
a lm ■ 



i=4 



Im 



(1) 



where a\ m is spherical harmonic transform of map and B\ is the 
beam function for the channel i supplied by the WMAP team. 
The condition that the CMB signal remains untouched during 
cleaning is encoded as the constraint [Wi][e] = [e] T [W] T = 1 
where [e] is a 4 x 1 column vector with unit elements. 

Following Tegmark & Efstathiou (1996), Tegmark et al. 
(2003) and Tegmark et al. (2000a) we obtain the opti- 
mum weights which combine 4 different frequency chan- 
nels subject to the constraint that CMB is untouched, [Wi] = 
[e] T [Cir7 ([e] T [Ci]- 1 [e]). Here the matrix [Q] is 



. m=t i j* 



(2) 



In practice, we smooth all the elements of the C/ using a 
moving average window over Al = 1 1 before deconvolving by 
the beam function. This avoids the possibility of an ocassional 
singular C; matrix. The entire cleaning procedure is automated 



and takes approximately 3 hours on a 16 alpha processor ma- 
chine to get the 48 cleaned maps. One of the cleaned maps, 
C8, is shown in the Fig. 1. In all the 48 maps some residual 
foreground contamination is visibly present along a small nar- 
row strip on the galactic plane. In a future publication, we as- 
sess the quality of foreground cleaning in these maps using the 
Bipolar power spectrum method (Hajian & Souradeep 2003 ; 
Hajian et al. 2004a; Hajian & Souradeep 2005) and compare 
them to others such as the internal Linear combination map 
(ILC) of WMAP. For the angular power estimation that follows, 
the Kp2 mask employed suffices to mask the contaminated re- 
gion in all the 48 maps. 




FIG. 1 . — The cleaned map C8 for the K, Q2, V2, W13 combination (units 
mK). Residual foreground contamination visible along the galactic plane are 
well within the Kp2 mask that is applied before extracting the angular power 
spectrum. 



2.2. Power Spectrum Estimation 

We obtain cross correlated power spectrum from these 
cleaned maps after applying a Kp2 mask. In choosing pairs 
Ci & Cj to be cross correlated, we ensure that no DA is com- 
mon between them. Figure 2 lists and plots the 24 cross 
power spectra for which the noise bias is zero. The cross 
power spectra are corrected for effect of the mask, the beam 
and pixel window. These are accounted for by de-biasing 
the pseudo-C/ estimate using the coupling (bias) matrix corre- 
sponding to the Kp2 mask and appropriate circularized beam 
transform (Hivon et al. 2002). Figure 2 plots the 24 cross 
power spectra (binned) individually. The spectra closely match 
each other for I < 540. The 24 cross power spectra are then 
combined with equal weights into a single 'Uniform aver- 
age' power spectrum 3 . We also estimate the residual power 
contamination in the 'Uniform average' power spectrum from 
the unresolved point sources on running through our anal- 
ysis the same source model used by WMAP team to cor- 
rect for this contaminant (Hinshaw et al. 2003a). The model 
is derived entirely within the WMAP data based on fluxes 
and spectra of 208 resolved point sources identified in the 
maps (Bennett et al. 2003b). The residual power from un- 
resolved point sources is a constant offset of ~ 140/iK 2 for 
/ > 400 (and negligible at lower /). This residual is much less 
than actual point source contamination in Q, KA or K band and 
intermediate between V and W band point source contamina- 
tion. It is noteworthy that the method significantly tempers the 
point source residual at large I that otherwise is oc I 2 in each 
map. The final power spectrum is binned in the same manner 



3 There exists the additional freedom to choose optimal weights for combining the 24 cross-power spectra which we do not discuss in this work. 



Angular power spectrum of WMAP CMB anisotropy 
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as the WMAP's published result for ease of comparison. 

2.3. Error Estimate on the Power Spectrum 

The errors on the final power spectrum are computed from 
110 random Monte Carlo simulations of CMB maps for ev- 
ery DA each with a realization based on corresponding WMAP 
noise map (available at LAMBDA) and diffuse foreground con- 
tamination. The common CMB signal in all the maps was 
based on a realization of the WMAP 'power law' best fit A- 
CDM model (Spergel et al. 2003). We use the publicly avail- 
able Planck Sky Model to simulate the contamination from the 
diffuse galactic (synchrotron, thermal dust and free-free) emis- 
sion at the WMAP frequencies. The CMB maps were smoothed 
by the beam function appropriate for each WMAP's detector. 
The set of DA maps corresponding to each realizations were 
passed through the same pipeline used for the real data. Aver- 
aging over the 110 power spectra we recover the model power 
spectrum, but for a hint of bias towards lower values in the low 
/ moments. For I = 2 and I = 3 the bias is -27.4% and -13.8% 
respectively. However this bias become negligible at higher I, 
e.g. at / = 22, it is only -0.8%. 

The standard deviation obtained from the diagonal elements 
of the co variance matrix is used as the error bars on the C/'s 
obtained from the data. (The beam uncertainty is not included 
here, but is deferred to future work where we also incorporate 
non-circular beam corrections (Mitra et al. 2004).) 



24 Cross-Power Spectra 
C1 ®CA12 C13®CA24 
C2SCA11 C14 0CA23 
C3SCA10 C15 0CA22 
C4®CA9 C16 0CA21 
C5 ® CA8 
C6 ® CA7 
C7 ® CA6 
C8 ® CA5 
C9 ® CA4 
C10 8CA3 




Multipole, I 

FIG. 2. — The 24 individual cross power spectra corresponding to the cross 
correlations listed in this figure are plotted. The cross spectra show very small 
dispersion for / <^540. The average power spectrum is plotted in red line and 
blue error bars. The multipole range is on log scale for / < 100, and linear, 
thereafter. 



3. RESULTS 

We obtain a 'Uniform average' power spectrum of the 24 
cross spectrum following the method mentioned in section 2.2. 
The black curve in Fig. 3 compares our results with the WMAP 
published power spectrum plotted with red line. The power 
spectra from two independent analysis are reasonably close. In 
case of 'Uniform average' a maximum difference of 92 /j,K 2 is 
observed only for octopole. For the large multipole range the 
difference is small and for I = 752 it is approximately 48 fiK 2 . 
This is well within the la error bar (510 fiK 2 ) obtained from 
the simulation. The small, but systematic excess, at large multi- 
poles is precisely resolved when our 'Uniform average' is cor- 
rected for the expected residual power from unresolved point 



source contamination described in §2.2. The point source cor- 
rected power spectrum is shown in black line in this figure. The 
difference of this power spectrum with the published WMAP 
estimate is shown in the bottom panel of Fig. 3. The differences 
are well within the la error-bars estimated from the simulations 
described in §2.3. 




Multipole, I 

FIG. 3. — The 'Uniform average power spectrum is plotted in green line. The 
black line shows residual unresolved point source corrected power spectrum. 
The error-bars are computed from the diagonal elements of the covariance 
matrix obtained from our simulation pipeline. Beyond / = 400 all the error 
bars are shifted by 81 = 10 to visually distinguish between error bars obtained 
from our method and the WMAP's published error bars. The published binned 
WMAP power spectrum is plotted in red line with error bars. The multipole 
range 2 < I < 100 are plotted in the log scale to show the small / behavior of 
the power spectrum. 

We find a suppression of power in the quadrupole and 
octopole moments consistent with WMAP published result. 
However, our quadrupole moment (146/i/T 2 ) is a little larger 
than the WMAP's quadrupole moment (123/i/T 2 ) and Octopole 
(455 fj,K 2 ) is less than WMAP's result (611/ztf 2 ). The 'Uni- 
form average' power spectrum does not show the 'bite' like 
feature present in WMAP's power spectrum at the first acoustic 
peak reported by WMAP (Hinshaw et al. 2003a). We perform 
a quadratic fit to the peaks and troughs of the binned spectrum 
similar to WMAP analysis (Page et al. 2003). For the residual 
point source corrected ('Uniform average') power spectrum we 
obtain the first acoustic peak at / = 219.8 ±0.8 (220.8 ±0.8) 
with amplitude A7} = 74.1 ±0.3^ (74.4 ± 0.3/^T), the sec- 
ond acoustic peak at I = 544 ± 17 (545 ±17) with amplitude 
AT, = 48.3 ± IJlfiK (49.6 ± 1.2fjK) and the first trough at 
I = 419.2 ±5.6 (418.7 ±5.5) with amplitude AT, = 41.7 ± 1 
^(42.2 ±0.9^). 

As cross checks of the method, we have carried out analysis 
with other possible combinations of the DA maps. 

1. The WMAP team also provide foreground cleaned 
maps corresponding to Ql to W4 DA (LAMBDA). 
The Galactic foreground signal, consisting of syn- 
chrotron, free-free, and dust emission, was re- 
moved using the 3-band, 5-parameter template fitting 
method (Bennett et al. 2003b). We also include K and 
Ka band maps which are not foreground cleaned. The 
resulting power spectrum from our analysis matches 
closely to the 'Uniform average' power spectrum. 

2. Excluding the K and Ka band from our analysis we get a 
power spectrum close to the 'Uniform average' results. 
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Notably, we find a more prominent notch at / = 4 similar 
to WMAP's published results. 

This is a clear demonstration that the blind approach to fore- 
ground cleaning is comparable in efficiency to that from tem- 
plate fitting methods and certainly adequate for a reliable esti- 
mation of the angular power spectrum. 

4. CONCLUSION 

The rapid improvement in the sensitivity and resolution of 
the CMB experiments has posed increasingly stringent require- 
ments on the level of separation and removal of the foreground 
contaminants. Standard approaches to foreground removal, 
usually incorporate the extra information about the foregrounds 
available at other frequencies, the spatial structure and distri- 
bution in constructing a foreground template at the frequencies 
of the CMB measurements. These approach could be suscep- 
tible to uncertainties and inadequacies of modeling involved in 
extrapolating from the frequency of observation to CMB obser- 
vations. 

We carry out an estimation of the CMB power spectrum from 
the WMAP first year data that is independent of foreground 
model and evades these uncertainties. The novelty is to make 
clean maps from the difference assemblies and exploit the lack 
of noise correlation between the independent channels to elim- 
inate noise bias. This is the first demonstration that the angular 



power spectrum of CMB anisotropy can be reliably estimated 
with precision solely from the WMAP data (difference assembly 
maps) without recourse to any external data. 

The understanding of polarized foreground contamination in 
CMB polarization maps is rather scarce. Hence modeling un- 
certainties could dominate the systematics error budget of con- 
ventional foreground cleaning. The blind approach extended to 
estimating polarization spectra after cleaning CMB polarization 
maps could prove to be particularly advantageous. 
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